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Jet vetoes play an important role at the LHC in the search for the Higgs and ultimately in precise 
measurements of its properties. Many Higgs analyses divide the cross section into exclusive jet bins 
to maximize the sensitivity in different production and decay channels. For a given jet category, the 
veto on additional jets introduces sensitivity to soft and collinear emissions, which causes logarithms 
in the perturbative expansion that need to be resummed to obtain precise predictions. We study the 
higher-order resummation properties of several conceptually distinct kinematic variables that can 
be used to veto jets in hadronic collisions. We consider two inclusive variables, the scalar sum over 
Pt and beam thrust, and two corresponding exclusive variables based on jet algorithms, namely the 
largest pr and largest beam thrust of a jet. The inclusive variables can, in principle, be resummed 
to higher orders. We show that for the jet-based variables, there are dual effects due to clustering 
in the jet algorithm for both large and small jet radius R that make a complete resummation at or 
beyond next-to- leading logarithmic order (NLL) challenging. For 7? ~ 1, the clustering of soft and 
collinear emissions gives 0(1) contributions starting at next-to-next-to-leading logarithmic order 
that are not reproduced by an all-orders soft-collinear factorization formula and therefore are not 
automatically resummed by it. For R <C 1, clustering induces logarithms of 7? that contribute 
at NLL in the exponent of the cross section, which cannot be resummed with currently available 
methods. We explicitly compute the leading jet clustering effects at 0(a 2 a ) and comment on their 
numerical size. 



I. INTRODUCTION 

The LHC is at the brink of discovering (or ruling out) 
a Standard- Model-like Higgs boson [TJ [5] . Many Higgs 
and new-physics analyses divide the data into categories 
based on the number of hard jets in the final state. 
This "jet-binning" significantly enhances the experimen- 
tal sensitivity to various combinations of Higgs produc- 
tion and decay channels. The experimental distinction 
between different channels is of course an important re- 
quirement to determine the Higgs properties. Any cou- 
pling measurement requires precise and reliable theory 
predictions of the jet cross sections used in the experi- 
mental analysis. 

By vetoing jets, one restricts the phase space for addi- 
tional emissions. This makes the cross section sensitive to 
soft and collinear radiation, which induces Sudakov dou- 
ble logarithms of the jet-veto variable at each order in 
perturbation theory. At small enough cuts the logarith- 
mic corrections become large and dominate the pertur- 
bative series, which degrades the reliability of fixed-order 
perturbation theory. To obtain precise predictions and 
robust uncertainty estimates it is important to under- 
stand the structure of large logarithms and ideally resum 
them to all orders in perturbation theory. The form of 
the jet-veto logarithms and the precise structure of the 
logarithmic series depends on the details of how the jet 
veto is imposed. 

In this paper we consider four kinematic variables that 
are conceptually and theoretically distinct, and study the 
prospects for resumming the respective large logarithms 
when these (types of) variables are used to veto central 
jets. For definiteness we will discuss the process gg — ¥ 



H + jets in the following, but our analysis applies to the 
production of any color-singlet in conjunction with a veto 
on central jets. It also carries over to generic processes 
where one requires N signal jets in the final state and 
places a veto on additional jets, i.e. when measuring an 
exclusive /V-jet cross section. 

In principle, one can think of various kinematic vari- 
ables that can be employed to enforce a veto on hard 
emissions. The four variables we will discuss are summa- 
rized in Table |l] and are described next. They are classi- 
fied as either inclusive or exclusive (jet-based) variables 
and according to their sensitivity to either the virtuality 
or pt of emissions. 

The inclusive variables we consider are beam thrust [3] 
and the scalar sum of pr, 

T B = Y}PTm\e-^- Y K (1) 
m 

E T = ^\PT m \, (2) 

m 

where f>Tm and y m are the transverse momenta and ra- 
pidities with respect to the beam axis of all particles in 
the final state, but excluding the Higgs decay products, 
and Y is the rapidity of the Higgs boson. 1 For these 
inclusive variables one sums over all hadrons m in the fi- 
nal state, which at the perturbative level corresponds to 
constraining the sum over all emissions. As a result such 



For cases like H — > WW where Y cannot be measured directly, 
one can also consider the analog of Tb defined in the hadronic 
center-of-mass frame by setting Y = in Eq. (TTj) . 
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inclusive variables are theoretically the cleanest and best 
understood. The complete resummation of logarithms at 
small Tb and Et has been carried out to next-to-next- 
to-leading logarithmic order (NNLL) [H [5] and next-to- 
leading logarithmic order (NLL) [5] , respectively, and the 
extension to higher logarithmic orders poses no concep- 
tual difficulties. The inclusive nature of these variables 
makes it challenging to measure them in a hadron-collider 
environment. These experimental issues can be mitigated 
by summing over jets rather than hadrons, as is done in 
current event-shape measurements at the LHC [7J [5] , or 
by only summing over charged tracks, which can then be 
corrected for in the measurement. 

The difference between Tb and Et is in the differ- 
ent rapidity weighting of emissions. In the rest frame 
of the Higgs, where Y = 0, Tb is equivalent to the 
sum over the small light-cone component of momenta, 
7b = J2k(Ek — \p%\) ~ t/rriH, where t is the spacelike 
virtuality of the colliding hard partons after initial-state 
radiation. Hence, by measuring beam thrust Tb one is 
sensitive to the virtuality scale of emissions. In contrast, 
by measuring Et one is sensitive to the pt scale of emis- 
sions, and one might call it "beam broadening". This 
difference in sensitivity to virtuality vs. pt causes the 
logarithmic series at small Tb or Et to have very differ- 
ent structures. In the parton shower, this difference is 
analogous to the different Sudakov form factors for vir- 
tuality and pt ordered showers. 

The exclusive jet-based variables we consider are the 
largest beam-thrust or pt of a jet, given by 

Tj= max \p Tm \ e -^- Y \ , (3) 
mej(R) 

p Tj = max \p Tm \ , (4) 

mSj'(-R) 

where pTj and yj are the jet's transverse momentum and 
rapidity. They have the same sensitivity to the virtuality 
or pt scale of emissions as their inclusive counterparts Tb 
and Et in Eqs. ([I]) and ([2j. The exclusive variables are 
based on identifying jets j(R) of size R and considering 
the largest contribution from a jet. Pcrturbatively, this 
corresponds to constraining the maximum of all emissions 
at a typical "resolution scale" i?, which is in contrast to 
constraining the sum of emissions as in the inclusive vari- 
ables. We focus on jets using the fcx class of algorithms 
[5HT3"]. The jet-based variables are more straightforward 
to use experimentally, and a pt veto on jets using the 
anti-ZcT algorithm is the common choice in experiments. 
On the other hand, the variable's inherent dependence on 
the jet algorithm and the resolution scale R make them 
considerably less tractable theoretically. The resumma- 
tion for pTj was considered recently in Refs. [141 f!5] . 

We implement the veto on central jets by putting a 
constraint 

Tb < T cut , E T < P c T ut , 

Tj < T cut , PTj < P c T ut , (5) 





sensitive to 




virtuality 


Pt 


inclusive 


Tb 


Et 


jet-based (exclusive) 


Tj 


PTj 


constraint 


^ ST- CUt 


S Pt 



TABLE I: Classification of different jet-veto variables. 

and define the small parameter A by 

A 2 = T cut /m H or A = p^ fm H ■ (6) 

The cross section with a T cut or p™* veto contains Su- 
dakov logarithms a™ ln m A with m < 2n. In the region 
of small A the logarithms dominate the cross section and 
we want to resum them. At the same time, any nonlog- 
arithmic contributions that depend on T cut or p™* are 
suppressed by relative powers of A, and can be added to 
the resummed result. 

For the jet-based variables, one can distinguish two 
cases of how the jet size R is counted relative to A, we 
can consider either i? ~ A or > A. 2 For this dis- 
tinction it is irrelevant whether one considers R ~ 1 or 
fl< 1. We stress that which of these two formal cases is 
in the end more appropriate in practice, for given numer- 
ical values of R and T cut / m -W or p^/niH, is a separate 
question that needs to be studied numerically and will 
be addressed in Sec. [V] 

As we will show in this paper, there are competing ef- 
fects in either limit arising from the dependence on the 
jet algorithm, which spoil the complete logarithmic re- 
summation. Jet algorithm dependent effects first arise 
at 0{a 2 s ) in the form 

a 2 /g(i?)mA. (7) 

The function /^(R) depends on the algorithm and the 
jet- veto variable, and in general contains terms of 0(R 2 ), 
constant terms, and lni? terms. These were calculated 
in Ref. [14] for the pTj veto. 

For R ^> A the jet algorithm mixes soft and collinear 
contributions in the measurement, giving rise to a term 
in Eq. ^ of the form 

a 2 i? 2 lnA, (8) 

at leading order in R. When keeping the full R depen- 
dence, these soft-collinear mixing terms inhibit an all- 
order soft-collinear factorization in the measurement at 



2 This counting is a natural proxy for the two general cases of 
R either scaling as a positive power of A or i? not scaling with 
A. The difference amounts to whether we formally count O(R) 
terms as power corrections in A or not and logarithms of R as 
In A or not. 
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leading power in A. In Ref. [T5], a factorization formula 
for the pTj veto was presented in the limit A <C R ~ 1. 
Their derivation does not account for the effect of these 
mixing terms to all orders, effectively assuming that they 
are power-suppressed, and hence breaks down for R ~ 1. 

For R ~ A, the mixing terms inhibiting the soft- 
collinear factorization of the measurement can be re- 
garded as power corrections in A, so factorization goes 
through. However, in this limit clustering effects that 
change the boundary of the jet at each order in a s intro- 
duce corrections that depend on In R. These are referred 
to as clustering logarithms and were first pointed out in 
Ref. [TS] and were studied in Refs. [P7H2T] . The cluster- 

(2) 

ing logarithms in /^ lg (R) give a contribution of the form 



a 



I In R In A . 



(9) 



They are distinct from those previously studied since 
they are associated with collinear rather than soft di- 
vergences within each jet. The 0(a 2 ) term in Eq. (m) is 
the first in an all-orders series of terms of the form 



a^fW(R) In A, 



fZ\R)Dln k R, k<n-l. (10) 



For R ~ A we have to count In R ~ In A, so these cluster- 
ing logarithms give a new contribution at NLL at each 
higher order in a s . Hence, they spoil the complete loga- 
rithmic resummation at NLL and beyond. 

At a formal level, this means one is stuck between a 
rock and a hard place. We would like to consider R 
as small to justify not resumming soft-collinear mixing 
terms, by treating them as unresummed power correc- 
tions. But at the same time we would like to treat R 
as large, to avoid having to count mi? terms as large 
logarithms. 

A possible way forward would be to resum those log- 
arithms whose all-order series is known while simulta- 
neously quantifying the effect of the terms that cannot 
be resummed. If these effects can be appropriately folded 
into uncertainty estimates, then a reliable theoretical pre- 
diction can be obtained. The formal power counting is 
less important in this case, and effectively the two formal 
limits for R are unified by reflecting the important terms 
from each case in the perturbative uncertainties. 

The remainder of this paper is organized as follows. In 
Sec. |TT] we present an overview of the necessary steps to 
achieve soft-collinear factorization and resummation, fo- 
cusing on the properties of the measurement. We present 
the factorization formulas for the inclusive observables in 
Sec. Ill Al and for the exclusive observables in Sec. Ill Bl In 



Sec. Ill we discuss the soft-collinear mixing in detail and 
show that it gives an 0(R 2 ) contribution to the rate. We 
demonstrate this by power counting as well as explicit 
calculation of the mixing terms at 0(a 2 ), the details of 
which are given in Appendix [B| In Sec. |IV| we discuss the 
clustering logarithms in detail and calculate their contri- 
bution at 0(a 2 ), with the details of the calculation given 
in Appendix [C] We use non-Abelian exponentiation and 



the collinear sensitivity of the clustering logarithms to 
show that their contribution to the cross section expo- 
nentiates and for R ~ A contributes at NLL in the expo- 
nent. In Sec. |V]we summarize our findings and give an 
outlook based on a numerical analysis of the size of the 
0{a 2 ) mixing and clustering terms for different values of 
R. 



II. FACTORIZATION FORMULAS 

The large Sudakov logarithms in the jet-vetoed cross 
sections arise as a remnant of the cancellation of soft and 
collinear IR divergences between virtual and real contri- 
butions. Their resummation is thus intimately tied to 
the universal structure of QCD amplitudes in the soft 
and collinear limit. Furthermore, resummation relies on 
the fact that soft gluon emissions from energetic particles 
are eikonal and that the total soft eikonal matrix element 
factorizes from the remaining amplitude. 

A convenient framework to study the logarithmic 
structure of the cross section is provided by soft-collinear 
effective theory (SCET) [2"5H27] . which makes the soft- 
collinear limit of QCD manifest at a Lagrangian and 
operator level using a systematic power expansion in 
the small parameter A. The resummation of Sudakov 
logarithms is then achieved by standard effective-theory 
methods through a systematic scale separation and renor- 
malization group evolution between the scales. For a de- 
tailed discussion of this procedure in the context of SCET 
we refer the reader to the literature. In the following, we 
give a schematic overview of the basic steps, concentrat- 
ing on the features that are important for our further 
analysis. 

After matching full QCD onto SCET, the cross section 
in SCET for gg — > H with no hard jets in the final states 
has the schematic form (for details see e.g. Refs. [21 [5]) 

dCTgg^H ~ \CggH{^)\ 2 {PaPb\Ogg H {lJ^ MOgg H {lj)\PaPb) ■ 

(11) 

The incoming (anti)protons have momenta 



Lt 

Pa,b 



E r 



with 



(1,5), 



(1,- 



(12) 

The Wilson coefficient C gg #(/i) arises from matching 
onto the O gg H{n) operator in SCET and encodes the 
hard-scattering contributions, which live at the hard 
scale fin ~ "iff- The measurement operator M imple- 
ments the phase-space cuts and measurements on the fi- 
nal state. The QCD dynamics in the soft and collinear 
limits, at leading order in the expansion parameter A, are 
encoded in the SCET operator matrix element. At this 
point, the renormalization group evolution (RGE) of the 
hard Wilson coefficient can be used to sum logarithms 
of the form ln(/i cs /m^) ~ In A that arise as the ratio of 
some low soft-collinear scale fi sc ~ Am# (set by the mea- 
surement in the matrix element) and the hard interaction 
scale ~ ma ■ 
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The matrix element in Eq. (Ill contains further log- 



arithms of A due to the different scaling of soft and 
collinear contributions. These logarithms are part of the 
full logarithmic structure of the cross section, and their 
resummation requires the separation of soft and collinear 
contributions. In SCET, this separation proceeds in two 
steps. First, the decoupling of soft emissions happens 
via a held redefinition of the collinear quark and gluon 
fields [25 , after which O gg n takes the form 

O ggH =HO a O s O b = H B» a± T [y\ a y nb ] B nb± M . (13) 

Here, H is the Higgs field, O a and Ob are collinear gluon 
fields (B n ±) in the forward (n a ) and backward (rib) beam 
directions, and O s is an operator of soft lightlike Wilson 
lines (y n ) along the n a and n& directions. After the soft- 
collinear decoupling, the SCET Lagrangian contains no 
interactions between soft and collinear fields at leading 
order in A. ____ 

Second, the measurement M must be separated into 
collinear and soft components that act independently on 
the soft and collinear final states. Schematically, 



M = M a x M b x Ms + SM. 



(14) 



The operators M a , Mb, and M s are obtained by re- 
stricting M to act only on rt Q -coilinear, rib-collinear, and 
soft fields respectively. The total contribution from the 
remainder 5M must be power suppressed in A (such that 
to all orders it can contribute at most terms that scale 
like A In" A). 

With these two ingredients, the matrix element in 
Eq. (11) factorizes into independent soft and collinear 



matrix elements, 



Ba{p) ~ ( Pa \OtM a O a \p a ){n) 

B b (fjL) ~ (pb\OtM b O b \p b )(ji), 
S(n)~(0\OiM a O a \0)fr). 



(15) 



Here S is a soft function and B is a beam function [3J [251 
|2T)] . which describes the collinear initial-state radiation 
from the parton entering the hard interaction, and can be 
calculated perturbatively by matching onto the standard 
parton distribution functions [see Eq. (B3)]. It follows 
that the cross section factorizes as well, 



da 



gg^H 



\C ggH (fi)\ 2 [B*(n) x fl 6 (/i) x S(n)] . (16) 



The collinear and soft matrix elements are renormal- 
ized objects and depend on a renormalization (or sep- 
aration) scale ijl. The RGE running between the natu- 
ral collinear and soft scales then resums logarithms of 
the form ln(/i c //z s ) ~ In A that are present in the soft- 
collinear matrix element itself. 

From this discussion, we can see that resumming the 
logarithmic series in In A requires an explicit separation 
of the degrees of freedom whose associated scale depends 
on A. In the effective-field theory context the resumma- 
tion is then achieved by RGE methods. The complete 



all-order resummation requires this separation to hold to 
all orders in perturbation theory. This is guaranteed for 
the hard-scattering factorization in Eq. (11), which es- 



sentially amounts to expanding QCD in the soft-collinear 
limit. It also holds for the soft-collinear operator decou- 
pling, which is independent of the measurement. How- 
ever, for the measurement function it means that the 



separation in Eq. ( 14 ) has to hold for any number of soft 
and collinear particles in the final state. As we will see 
below, this requirement is satisfied by the inclusive ob- 
servables, but for the jet-based observables it provides a 
nontrivial constraint and only holds for R ~ A. 

In the following we discuss the soft-collinear factoriza- 
tion properties of our four observables and the resulting 
factorization formulas for the cross section. We consider 
the inclusive observables in Sec . |II A| and the jet-based 
ones in Sec. |IIB| In Appendix |A| we discuss the RGE 
constraints on the factorization formulas for these ob- 
servables. 



Inclusive Observables 



The full measurement operator M entering Eq. (11) 
directly follows from the definitions of the observables in 
Eqs. (JlJ and (|2). For example, for beam thrust its action 
on a given fina state with a set of momenta {p m } is 

M(T B ) = s(Tb - J2\PTm\ e-l» — ^1) , (17) 



and analogously for Et- Since the inclusive observ- 
ables simply sum over particles, we can write them as 
a sum over separate contributions from n a -collinear, n b - 
collinear, and soft particles, 



Tb — Tsa + Tsb + Tbs , 
Et = Eto + Exb + Ets ■ 



(18) 



This means there is no contribution SM at leading power 
which mixes the different sectors, and we can factorize 
the measurement function in terms of a convolution, 

M(k) = J dk a dk b dk s S(k — k a — k b — k s ) 

x M a (k a ) M b {kb) M s {k s ) , (19) 

where k here stands for either Tb or Et- In the case of 
Tb we have 



M a (TBa) = S(r Ba - J2 \PTm\e Y - V " 



m£n a — coll 



Vm-Y 



M a (T B b) = ^{TBb - \PTm\e 
m(E.rib — coll 

M s {T Bs ) = 6(Tbs - Yl e" l!/m - y| ) • (20) 

rn^soft 
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Note that for the collinear sectors, the contribution from 
the respective opposite hemisphere does not contribute, 
because it is power-suppressed. Similarly, for Et we have 



M a {E Ta ) = 5(E Ta - 
M a {E Tb )=5(E Tb - 
M s (E Ts )=5[e Ts - 



E 

m£n a — coll 

E 

m£nb— coll 



\PTm 
\PTm 



E 

rn£ soft 



(21) 



Using n a and n b we define light-cone coordinates, 



{n a -p,n b -p,p ± ) 



P^ = n .p|+n t .p^ 



-p1- 

(22) 

In terms of these, n a -collinear particles have momentum 
scaling p a ~ m#(A 2 , 1, A), while nt>-collinear particles 
have momentum scaling p;, ~ mjj(l, A 2 , A). The soft and 
collinear contributions to Tb can be written as 



Tbs = 



e n. 



Pn 



Ti 



Bb 



E 

m£soft 



/ min{e 



e Y n b ■ P b . 
Y n b -p m } . 



(23) 



where P a>b is the total momentum of all n a ^-collinear 
final-state particles. Since Tb measures the small light- 
cone components, soft particles that contribute to the 
measurement of Tb have ultrasoft (us) momentum scal- 
ing p us ~ itih (A 2 , A 2 , A 2 ). The appropriate version of 
SCET for this case is called SCETi. Since p 2 ab ~ 
\ 2 m 2 H 3> p 2 s ~ A 4 m 2 ? , soft and collinear degrees of free- 
dom in SCETi are separated in virtuality and the RGE 
running in SCETi describes evolution in invariant mass. 

The factorization formula for the beam thrust distri- 
bution is [3j [5] 



da 



a H gg (m H ,p,) JdY J dk a dk b 

x B g (m H k a ,x ai p) Bg(m, H k b ,x b ,fi) 
x S b 3 {Tb ~k a - kb,fi) , 



where 



and 



m H Y 

— — e . 



x b 



0-Q = 



V2G 



(24) 



(25) 



(26) 



5767r£ c 2 n ■ 

The convolution between beam and soft functions in 



Eq. (24 1 is a direct consequence of the convolution of 



the measurement in Eq. ( 19 ) 



In contrast to Tb, Et measures the perpendicular com- 
ponent of momentum, which means soft particles con- 
tributing to the measurement of Et have soft momen- 
tum scaling p s ~ m#(A, A, A). The appropriate version 



of SCET for this case is called SCET„. In SCET n , soft 
and collinear degrees of freedom have the same virtual- 



ity scaling, p 2 b ~ pj 



X 2 



H' 



but are still separated 



in rapidity. As a result, the RGE running in SCETn 
describes both evolution in invariant mass as well as ra- 
pidity [301 EI]. 

The factorization formula for the Et distribution is [6] 



da 
dE T 



a Hgg(m Hl fi) JdY J dk a dk b 

x B g (m H , k a ,x a , fi, v) B g (m H ,k b , x b , (j,, v) 



x S?f \E T — k a — k b , (j,, v 



(27) 



The soft and beam functions are different from those 



in Eq. ( 24 1 because they contain a different measure- 



ment. Beam functions with p^-dependence were studied 
in SCET in Refs. [15l EDH36] . Since the beam and soft 



functions in Eq. (27 1 are renormalized and RG evolved 



in rapidity, they depend on another scale v from this 
running. 

Higher-order resummation at small Tb and Et can be 
carried out systematically using Eqs. (24) and (27). Ex- 
tending the Tb resummation to N 3 LL requires determin- 
ing the two-loop beam and soft functions, while resuming 
Et to NNLL requires determining the two-loop beam and 
soft non-cusp anomalous dimensions. 



B. Exclusive Observables 

Vetoes on individual jets rather than an inclusive veto 
on the final state requires a more careful understanding 
of the role of the jet algorithm. When putting a cut on Tj 
or p Tj in Eqs. ^ and Q, the full measurement operator 
acting on the complete final state in Eq. (11 ) is 



x jGt (T cut ) = [] e (\p Tm \e-\ y ™- Y \ < T cut ^ 

mej(R) 

M et (p T at )= LI 9 (\PTm\ <PT Ut ) , (28) 

where we consider jets j(R) defined by the k^ class of 
algorithms. Note that integrating Tj or pTj up to a 
cut turns the maximum condition in their definitions in 
Eqs. (|3) and Q into the simple product of 9 functions 
in Eq. (281). Hence, for the exclusive variables it is more 
convenient to consider the integrated cross section with 
a cut rather than the differential spectrum as in the in- 
clusive variables. 

We now want to separate the measurement function 
into components that act independently on soft and 
collinear final states. That is, we want to express the 
full jet veto in Eq. |28|) as 



M iet (k c 



: A4 J a et (fc cut )A4 J 6 et (/c c 
+ SM' 3Ct (k cut ), 



jet / 7„CUt ~i 



)Mf(k 



(29) 



G 



where here and in the remainder of this section, k cut 
stands for either T cut or p^ . The functions Mf\k cut ) 
for i — a,b, s arc defined by the full measurement applied 
to n a -collinear, nf,-collinear, and soft particles, respec- 
tively, which also defines the remainder 5AA lct (k cut ). For 
this separation to be meaningful, the remainder should be 
power-suppressed, which requires that the full measure- 
ment M 3et does not mix constraints on collinear and soft 
particles. This is a nontrivial condition, since the veto 
on any individual jet is not allowed to mix constraints 
between sectors. If a jet has a collinear component p c 
and a soft component p s , then the veto condition, e.g. 
for p™' , is 



Since M sct in Eq. (29) is a simple product, and not a 



\PTc+PTs\ <P°T t 



(30) 



This prevents M lct from factorizing into separate soft 
and collinear components and gives a contribution to 
8M^ ct . Therefore, to preserve factorization this scenario 
should only happen with a power-suppressed rate, such 
that at leading power each jet contains either only soft or 
only collinear final states. In that case, we can perform 
the veto separately on jets in each sector, i.e., we can 
write the product over all jets in Eq. (28 1 as products 
over soft and collinear jets as in Eq. (29). 

We shall show that for R ~ A, 5M iCt indeed gives a 
power-suppressed contribution to the rate. In this limit, 
the following factorization formula holds for the Tj veto: 

<r(T cut ) =<T H gg (m H ,v) J dYBf (m H T cut ,w) 

x B'f (m H T cu \x bl M ) Sf(T cu \ M) , (31) 



where x a ,Xb and <7o a rc defined in Eqs. (25) and (26). In 
the same limit, an analogous factorization formula holds 
for the pxj veto: 

°(PT Ut ) = <ToH gg {m H iH) J dYBf (m H ,p c ^,x a , fi,v) 



B'f (m H , p§?* , £&, fi, v) Sf g (p§P* , M, v 



(32) 



An equivalent form of this factorization formula was de- 
rived in Ref. [TS] , and the NLL resummation for p™* was 
performed in Refs. fHl ITS] . 

In Ref. [TJ] jet-algorithm dependent effects were cal- 
culated at fixed 0{a 2 s ). These results were interpreted 
in Ref. [15] in terms of a two-loop anomalous dimension 
and used to extend the resummation based on Eq. ( 32 ) 
to NNLL, working in the limit A <C R ~ 1 to avoid In R 
clustering logarithms. However, we will show explicitly 
in the following that for R ~ 1 Eqs. (31 1 and (32) do not 



reproduce the all-order structure of QCD beyond NLL. 
We will also see that for R ~ A, where the factorization 
formula holds, only some parts of the 0(a 2 s ) contribu- 
tions from Ref. j!4j are correctly interpreted in terms 
of anomalous dimensions and used in conjunction with 
Eq. pi 



convolution as for the inclusive variables in Eq. (19), the 



factorized cross sections in Eqs. (31) and (32) now con 



tain a product of beam and soft functions rather than 
a convolution. Also, each function explicitly depends on 
the jet algorithm used, in addition to the jet-veto vari- 
able itself, and includes independent clustering effects. 
Note that the 0(a s ) results do not yet depend on the 
effects of the jet algorithm. Nevertheless, the resummed 
cross sections are different for the inclusive and exclu- 
sive observables starting at NLL because of the different 
structures of their factorization theorems. This reflects 
the fact that constraining the sum of emissions provides 
a very different phase-space constraint than constraining 
each individual emission for more than one emission. 

To understand the role of the jet algorithm in vetoing 
on individual jets and how it impacts the logarithmic 
series, it is useful to express the measurement function 
in the form 

M ict = {M a + AMf) (M b + AMf) {M s + AMf) 
+ 6M ict , (33) 

where AMf is defined to contain the jet-algorithm de- 
pendence within each of the collinear and soft sectors, 



Mf\k c 



M t {k c 



AMf\k c 



(34) 



for i = a,b, s. The definition of AMf is subtle, since it 
depends on what we define the corrections due to cluster- 
ing relative to; namely it depends on the precise choice 
of Mi, which is independent of R and the jet algorithm. 
To study the effect of clustering we choose Mi to be the 
inclusive Tb or Et measurement, 



M t (k c 



E 



m£i tn —sector 



(35) 



where k m is the Tb or Ex contribution from each particle. 

At 0(a s ) there are no jet-algorithm effects, since there 
is only a single, either soft or collinear, final-state par- 
ticle. This means that ^A^^ 1 -* and AM^ are zero, and 
M* et ^' reduces to a sum over terms with one nontrivial 
constraint in each sector, 



M jct{1) (k c 



6(ki < k cut ) 



(36) 



This gives the same 0(a s ) contribution as the integral 
over the k = T or Et distribution with k < k cut . 

Starting at 0(a 2 ) the role of the jet algorithm must be 
understood. When R A, soft-collinear mixing effects 
are important. They give a correction to the cross section 
of the form 



8af c - \C ggH {^)\ 2 

x (p a Pb\O ggH (^SM ic 



O ggH {^)\PaPb) , (37) 
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(a) 



(b) 



FIG. 1: (a) The contribution to soft-collinear mixing at ct 2 s 
from the independent emission of a collinear and soft gluon 
clustered into a single jet, relevant when R 3> A. (b) Cluster- 
ing corrections in the beam and soft functions from correlated 
emissions, relevant when R ~ A. 



where the measurement operator SA4 }Qt is defined by 
Eq. ( 29 ) as the difference between the full measurement 

At 



and its restriction to the soft and collinear sectors. 
O(al), the soft-collinear mixing contribution arises from 
the independent emission of a collinear and soft particle 
as shown in Fig. [lja). It is given by the measurement 
function 

SM & (2) ( fc cut) = 6M jct (2) (fc cu t) + SM ict (2) (fcCut) _ ^ 

where the two terms correspond to a soft particle clus- 
tering with either a n a -collinear or n;,-colhnear particle. 
For the fcx-class algorithms, which cluster particles with 
AR < R, 



5M% {2) {k mt ) = 9(AR as < R) [9(k ict < fc cut ) 
- 6(k a < k cut )d(k s < fc cut )] , 



(39) 



where fcj C t is the value of the jet observable after com- 
bining the soft and collinear particles. Since SA4^ ct is 
defined as the correction to the factorized measurement, 
the second term subtracts the corresponding 0(a 2 ) con- 
tribution from Ma*^ M^* - ^ ■ 



In Sec. Ill we calculate the (leading) contribution from 
Eq. ( 39 1 to the cross section explicitly, and we will see 
that it has the form ~ a 2 R 2 ln(m# /k cut ) for both 7} and 
PTj vetoes. There are also factorized clustering correc- 
tions of the same form in AAif 1 (for i = a,b, s) from 
two independent emissions within each sector, and which 
are part of the two-loop soft and beam functions. The 
total QCD contribution from independent emissions thus 
arises as the sum of the factorized contributions and the 
soft-collinear mixing contribution. This shows explicitly 
that the factorized cross sections in Eqs. (31 1 and (32) do 



not reproduce the full NNLL structure of QCD for gen- 
eral R ~ 1, since they do not contain the soft-collinear 
mixing contributions. 

When R ~ A, clustering corrections from indepen- 
dent emissions, including the soft-collinear mixing term 
SA4^ Qt , can be regarded as power corrections. In this limit 
a different type of clustering corrections in AMf' be- 
come important, arising from correlated emissions within 
the beam and soft functions, as shown in Fig. [ljb). At 
0{a 2 s ) they are described by the measurement function 



= y(AR u < R)6{k jct < k cnt ) 
+ 6(AR 12 > R)9(h < k cut )9{k 2 < fc cut ) 

- e(k t + k 2 < k cut ) , (40) 

where particles 1 and 2 both belong to sector i. In 



Sec. IV we calculate the contributions from Eq. (40), and 
we will see that they have the form a 2 In R\n(rriH~/ k cut ) . 
We also discuss the higher-order structure of clustering 
logarithms and argue that at 0(a") terms of the form 
a™ In 11-1 R\n(rnH /k cut ) contribute in the exponent of the 
cross section. For R ~ A these terms are NLL. Since the 
coefficients of the clustering logarithms at each order con- 
tain a genuinely new algorithm-dependent contribution, 
which will generically be unrelated to lower orders, they 
cannot be resummed with present methods [21] . 



1. Rapidity Cutoffs on Jets 

In an experiment, the physical limitations of the de- 
tector impose a rapidity cutoff on measured jets. The 
parton luminosities naturally suppress forward, high-p-r 
jets, but in experimental analyses even moderate rapid- 
ity cutoffs (e.g., a cutoff of y cu t = 2.5) are used. Monte 
Carlo studies of the dependence on the cutoff with a pr 
veto on jets find a negligible effect for y cut ~ 4.5 but an 
0(10%) effect for y cut ~ 2.5 in the range of typical p™* 
values [51 IT4"]. 

The rapidity cutoff y cu t regulates the rapidity diver- 
gences present in the soft and beam functions for the PTj 
veto. In the bare functions, this effectively converts ra- 
pidity divergences into factors of y cu t- At 0(a s ), the RG 
structure of the hard, beam, and soft functions implies 
that these divergences have no effect on the fixed-order 
logarithms. The first order that the effect of the rapidity 
cutoff on the resummation can be observed is at 0(a 2 ), 
and, in principle, the cutoff could affect the 0(a 2 In 2 A) 
and 0{a 2 In A) terms. Although we do not consider these 
effects here, it would be interesting to study their impact 
on the logarithmic structure and resummation in more 
detail. 



III. SOFT-COLLINEAR MIXING 

The jet algorithm used to define Tj and pxj gives rise 
to jets that can contain both soft 3 and collinear particles, 
and the collinear beam radiation can give rise to multiple 



AM 



J ct (2)^cut^ 



s Our discussion in this section is mostly insensitive to whether the 
soft radiation is described by ultrasoft modes (a 7} veto, SCETj) 
or soft modes (a pxj veto, SCETn), but we will point out when 
subtleties arise. 
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jets. Placing a veto on such jets mixes the phase space 
constraints in each sector through conditions of the form 



T c + Ts< T 



\PTc+PTs\ <PT 



(41) 



and this mixing contributes to 5Ai ]et in Eq. (29). In this 



section we show that these soft-collinear mixing contri- 
butions scale like R 2 , and so R must scale as A for these 
effects to be power suppressed. 

Clustering jet algorithms build jets by merging parti- 
cles in the final state according to a distance metric pij 
between particles and a metric pi for single particles. In 
each clustering step, the minimum metric determines the 
action of the algorithm. If the minimum is a pairwise 
metric then that pair is merged into a new particle, and 
if the minimum is a single particle then that particle is 
promoted to a jet. For the anti-fcx algorithm [T3], these 
metrics are 



Pij = "fe ,Pt 



R 



(42) 



The algorithm will cluster particles with separation 
ARij < R together. 4 

Phase-space constraints from jet algorithms have been 
studied in SCET in Refs. [57PD] . By applying the canon- 
ical scaling of soft and collinear particles to the anti-fcT 
metrics in Eq. ( 42 1 we can determine the typical behavior 



of the algorithm. For example, for the case of a veto on 
the jet px, the metrics scale as 



R 



R A 



A" 



(43) 



This naive power counting suggests that clustering of 
soft and collinear particles will typically not occur, since 
soft particles are at central rapidities and collinear par- 
ticles are at large rapidities and as a result p sc 3> 
{Pea Pss, Pa Ps} (f° r ~ 1 or smaller). The ordering 
p sc > p ss , p cc implies that collinear and soft particles 
will separately cluster among themselves before cluster- 
ing with each other. The ordering p sc > p Cl p s implies 
that (groups of) collinear and soft particles will be pro- 
moted to jets before any soft-collinear clustering between 
them can take place. A similar argument applies to mea- 
suring 7j. 

Power counting along these lines, counting A < i? < 
ln(l/A), was used in Ref. [15] to argue that, at leading 
power, soft and collinear particles do not occupy the same 
jet. However, the canonical scaling in Eq. ( [43] ) is not suf- 
ficient to demonstrate that soft collinear clustering gives 
only a power-suppressed contribution to the cross sec- 
tion [41] . In particular, there is a collinear enhancement 



4 In some cases, a clustering sequence can pull two particles farther 
than R away from each other. Such configurations tend to be 



ooooo 




FIG. 2: The shape of jets in two spaces: y ~ </> (left), the 
geometry where jets are found, and 9 — (f> (right), which is 
the physical geometry. As jets become forward, their size in 
6 shrinks. In the right figure, the dashed blue line shows the 
approximate shape of jets if ARij is replaced with the physical 
angle 8ij between particles. 



in soft emissions along the beam directions (i.e., the di- 
rection of the soft Wilson lines). This is straightforward 
to see from the matrix element for eikonal emission, 



Ankonal ~ dy d(f> . 

Pt 



(44) 



The matrix element is flat in rapidity, meaning soft parti- 
cles will populate jets at a forward rapidities at an equal 
rate to jets at central rapidities. For a jet j in a given di- 
rection, the rate for a soft particle to be radiated within 
the cone of the jet scales as 



9(AR ja < i?) - R 2 



(45) 



If R ~ A, this rate is power suppressed and can be ne- 
glected. We will see this scaling confirmed in the explicit 
calculations of soft-collinear mixing below. 

The geometry of jets at a hadron collider changes sig- 
nificantly in the forward region. For central rapidities, 
the clustering condition ARij < R is approximated by 
Oij < i?, where 8 is the physical angle between i and j. 
At forward rapidities, the jet subtends a smaller angle Rg 
with respect to the beam axis. The size of the jets in <fi 
remains unchanged, so jets become elongated at forward 
rapidities. A comparison between jets in y — (f> space and 
9 — 4> space is shown in Fig. [2j If a jet is at rapidity yj 
with a corresponding angle to the beam Oj, then 



Rg w 26 j sinh R . 



(46) 



uncommon. 



At the characteristic angle 6 C ~ A of collinear emissions, 
the size of the jet scales approximately as XR. This means 
that collinear initial-state radiation from the incoming 
hard partons can create multiple jets of size R in the 
final state, since the size of the jet is smaller than (po- 
tentially parametrically smaller than) the total size of 
the collinear beam sector, which is an angle of order A 
around the beam direction. This is important because 
if each collinear sector created a single jet (around the 
beam), a veto on jets containing both soft and collinear 
particles would require SAi^ ct to only depend on the total 
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collinear momentum and the measurement could there- 
fore be expressed in a factorizable way. However this is 
not the case here. From this analysis, we see that jets 
can contain both soft and collinear particles. Vetoing on 
such jets requires a constraint on the combined soft and 
collinear contribution to the veto variable, as in Eq. (41 1. 



Therefore, soft-collinear factorization will be inhibited by 
these jets for R ~ 1. 

We can calculate the correction to the cross section 
from soft-collinear mixing, Sa^c in Eq. ([37]), for both 
jet-based veto observables at 0(a 2 ). The mixing cor- 
rection is given by the measurement function 5M) et ^ 



in Eqs. (38) and (39) inserted into matrix elements of 



OggH with one soft and one collinear emission. Because 
of the soft-collinear decoupling this is simply the prod- 
uct of the 0(a s ) matrix elements for one soft and one 
collinear emission. 

The form of the soft-collinear mixing terms are the 
same as the 0{a 2 s ) corrections when two independent 
emissions within either the soft or collinear sectors clus- 
ter (which is contained in AAif* 1 ). Therefore, the total 
correction from clustering of two independent emissions 
to the cross section is given by 



A(T indcp(2) = A(j mdcp(2) 



Act 



indcp (2) 
CC 



da 



jet (2) 
SC 



(47) 



where SS, CC, and SC denote the independent emission 
contributions from the soft sector, collinear sector, and 
from soft-collinear mixing, respectively. In addition to 
the soft-collinear mixing term, we also calculate the cor- 
rection from the soft sector, A<7oo ep . 5 The calculations 
are given in Appendix[Bj where we work to leading power 
in R, 0(R 2 ), and drop corrections of 0(R 4 ). 

The Tj veto is easier to calculate, as dimensional regu- 
larization fully handles all divergences. The soft-collinear 
mixing, Sa 3 ^, is the sum of contributions from each 
beam, So^q + Scr'g^ and is given in Eq. (B27). This 
includes a nontrivial zero-bin subtraction in the collinear 
sector. The two contributions Sa 3 gQ and Sct^q are iden- 
tical, and their sum is 



'SC 



>{T C 



(^) 



-clo 



-if In 



(48) 



7T y 3 T cut 
Here <7lo is the leading-order cross section, which is 
given in Eq. (B6). We can see explicitly that the soft- 



collinear mixing contains a single logarithm and scales 
as R 2 , as expected from the naive scaling calculation in 
Eq. ( 45 1 . The soft independent emission contribution, 



see Eq. (|B12j), is 

ACT mdep(2) (rcut) 



1 / aA y 
el 7r / 



V 



cut 2 



2 c 



12 



.-if. 
(49) 



The remaining contribution at this order is from the in- 
dependent emission of two collinear particles, Aer™p ep , 
whose form can be constrained by two facts. First, since 
the effective theory reproduces the IR structure of QCD, 
the collinear contribution must cancel the 1/e depen- 
dence in AfTg^ cp . Second, the scale dependence of the 
collinear matrix element can be determined by power 
counting. It follows that the A<7pp Cp contribution must 
have the form 



Act 



indcp (2) ^cut 
CC 



2e 2 

12' 



if. 

> 

'(50) 



Combining all three contributions, the leading in- 
dependent clustering effect from independent emissions 
with a cut on Tj is given by 



indcp (2) frj-cut\ 



(T cut ) = -cr LO 



°^)% R * ln pL (51) 
7T J 2 T cut 



For the pTj veto, the total 0{a 2 ) clustering effect from 
independent emissions in QCD, which includes the soft- 
collinear contributions, was calculated in Ref. |14j and 
found to be 



A(J indcp(2)^™t) = 



7T / 3 pSp* 



(52) 



This result will serve as a partial cross check on our re- 
sults. To calculate the soft-collinear mixing contribution 
for the pxj veto, we have to regulate rapidity divergences. 
We use the analytic regulator [42] in this case, for which 
v plays the same role as fi in dimensional regularization, 
and a the role of e. The soft-collinear mixing contribu- 
tions are different for Sa 3 ^ and Sa 3 gQ b , because of the 
asymmetry in the regulator between the two collinear 
sectors. We find 



jet (2) 
SC a 



CLCr 



a s C A ^ 



6 



-if 



2a 



T i / \m H J 
a s C A \2 2TT 2 



2a 



-if ln^ 

3 pip* 



(53) 



As for the Tj veto, the soft-collinear mixing for the pTj 
veto contains a single logarithm and scales as R 2 . The 
soft independent emission contribution for pTj is scaleless 
and thus vanishes 



A4 n s Cp(2) (PT Ut ) 



0. 



(54) 



5 Determining the O(a^) beam function matrix elements is sig- 
nificantly more involved, and is not required to demonstrate the 
presence of soft-collinear mixing. It is therefore beyond the scope 
of this work. 



Hence, to reproduce the full independent emission result 



in Eq. (52), the collinear contribution must be 

■a s C~^- 2 



Act 



indcp (2)/ cu t n 

CC \Pt I 



„cut 

ft 



(55) 
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The soft-collinear mixing contributions in Eqs. (48) 
and (531 are not isolated to 0(a 2 ). At higher orders, ad- 



ditional emissions can generate a tower of Sudakov dou- 
ble logarithms on top of these results and also generate 
higher-order mixing effects. Accurate uncertainty esti- 
mates require a better understanding of the soft-collinear 
mixing terms at all orders. 



IV. CLUSTERING LOGARITHMS 

We now turn to the form of clustering logarithms in 
the cross section. When i? ~ A, the cross section a(k cut ) 



satisfies the factorization formulas in Eqs. (31) and (32) 



Clustering of final state soft and collinear particles can 
give rise to logarithms of i? that become important when 
R <C 1. These logarithms arise as the remnant of a 
collinear divergence between particles in the jet, and are 
associated with connected webs (c-webs) in the matrix 
element [43]l44]. 

To see this, consider a two-parton final state with a 
PTj veto. The measurement function can be written as a 
constraint on each particle's px plus a correction factor 
for when they are clustered into a jet, 



AM 



jet (2) 



(56) 



where 



AM 



jet (2) 



e( P Tl <PT Ut )0( P T2 <PT Ut ), 
8(AR 12 < R) [9(\p T1 +p T 2\< Pt') 

-6(p T1 <p c T ut )6( PT2 <p c T ut )]. (57) 



Each individual measurement, M^J to and AM^to \ ^ s 
separately IR unsafe if the matrix element has a collinear 
singularity between the two particles. Since the total 
measurement, M^ ct ( 2 \ is IR safe, the divergence cancels. 
The remnant is logarithmic sensitivity to i?, schemati- 
cally 



a41 



i 



AMI 



- X -R- 



M {2) +AA/P t(2) 

■ /v 'vcto ' ZA - /vl veto 



lni?. 



(58) 



This structure persists for more particles in the final 
state. In general, if there are n final-state particles, there 
are at most n — 1 collinear singularities between them, 
each of which leads to a factor of In R from clustering ef- 
fects. We will show that the general form of the leading 
clustering logarithm at C(a") is 



c£>QnR) In 



„cut ' 

ft 



(59) 



where C^) (In R) contains at most n — 1 logarithms of R 
and the same form holds for T cut . 



A. O(ccg) Clustering Logarithms 



Although the division in Eq. ( 56 1 helps determine the 



order of the clustering logarithms, it is not well suited 
to define the clustering correction. Instead, we use the 



division in Sec. II B defining the clustering logarithms 
relative to the inclusive measurement. 

At 0(a 2 ), the clustering logarithms for the px veto 
have been calculated in Ref. [13]. We perform the soft 
function clustering calculation for both pTj and Tj in 
Appendix[C] since it is instructive to see the RG structure 
explicitly. The bare soft function corrections for Tj are 
given in Eq. ( |C18[ ) and for pTj in Eq. ( C17[ ), where for 
the latter we use the rapidity regulator [301131] . The bare 
corrections are UV-divergent and give a contribution to 
the two-loop soft anomalous dimensions of 

A^\T™\,)=(^) 2 C?(lnR), 

A 7 ^(^,,)=(^) 2 C( 2 T )(lni?), (60) 

where the coefficients are given by 

131 - 12tt 2 - 132 In 2 



4 2) (lni?) = Ci(: 



18 



lni? -0.936 



+ C A T R n f (- 23 g 41n2 lni? + 0.748 
131 - 12tt 2 - 132 In 2 



C$(lnR) = C^(— — In /,' - 1.12 



+ 



CaTrti s f- 



18 

23 - 241n2 
9 



lni? + 0.764 



(61) 



For the px veto, the contribution to the cross section 
stemming from this coefficient agrees with the results in 
Ref. [14] for the lni? term. The constant terms depend 
the observable the clustering effect is defined relative to. 
We use the inclusive Ex measurement, as in Eq. (40), for 



this purpose, whereas Ref. [TJ] uses the total px, which 
for two particles is px = \pti + PT2\- Either choice is 
possible and we have checked that we reproduce the con- 
stant terms in Ref. [14] when alternatively using their 
definition of the clustering correction. 

The soft anomalous dimension from these clustering 
effects must be canceled by the anomalous dimensions 
of the beam functions, since the total O gg H matrix ele- 
ment has only UV divergences that match the hard func- 
tion, and which are unrelated to clustering. It is worth- 
while to note that since no collinear singularities exist 
between soft and collinear particles, the cancellation of 
the divergences in the soft function from clustering must 
come entirely from the beam functions. It cannot come 
from a soft-collinear mixing term, which is power sup- 
pressed when i? ~ A. The cancellation of the beam and 
soft anomalous dimensions from clustering leaves a fixed- 
order contribution that has a logarithm of the ratio of the 
beam and soft scales. Note that the scale accompanying 
fi or v in the fixed-order beam and soft functions is fixed 
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by power counting. For the px veto, the beam v scale 
is rriH and the soft v scale is p™\ while for the T veto, 
the beam scale is \JmnT mt and the soft scale is T cut . 
Hence, the clustering effect at 0(a 2 ) in the cross section 
is 



A(J ( 2 )(T cut ) 

Aa (2) (PT Ut ) 



0LO 



V 7T 



ln"^C( 2 )(lni?). 

„cut Pt v I 
FT 



(62) 



This form matches that in Eq. (59). Because these terms 



are connected with the anomalous dimension, it is possi- 
ble to resum the logarithm of \n(m,H /k cut ) they contain. 
Since the anomalous dimension contributions from clus- 
tering involve only the beam and soft functions, the total 
evolution factor from RG evolution for each veto observ- 
able is 



#1 (Ast.Mj) 



U { l\v s , vj ) 



— exp 



exp 



^dM/o^)\ c ( 2) 

MS M V * J 



(63) 



Although the logarithms of R are exponentiated here, 
they are not connected with any scale in the effective 
theory. Below, we will show that the logarithms of R 
indeed appear in the exponent, meaning the above pro- 
cedure is correct. However, the non-cusp anomalous di- 
mension at O(ot™) contains a term ~ ln n .R, so for small 
R the perturbative series for the anomalous dimensions 
contains large unrcsummed logarithms. Equivalently, the 
perturbative series in the exponent is not resummed and 
contains terms a™ ln" _1 i?ln(mfl-/A: cut ). For R ~ A these 
are unresummed NLL corrections. Hence, the exponenti- 
ation of the 0(a 2 ) clustering logarithms is doing nothing 
to tame the NLL clustering logarithms at higher orders. 
The complete resummation of the clustering logarithms 
would require to calculate the entire tower of coefficients 
simultaneously. 



B. Higher-Order Structure of Clustering 
Logarithms 

The higher-order clustering logarithms are simplest to 
examine in the soft function, and RG constraints can be 
used to relate them to the beam functions and determine 
the overall effect on the cross section. The soft function 
at C(a™ ) can be written as 



S (n \ P T) = / &$> n A n {<S> n )M n {<5> n ,pT) 



(64) 



where <£>„ and A n are the n-particle phase space and 
matrix element. Non-Abelian exponentiation [43H46j im- 
plies that the eikonal matrix elements exponentiate, and 
can be factorized into c-web contributions, where a c-web 



is a diagram connecting eikonal lines that cannot be sep- 
arated into lower order c-webs by cutting each eikonal 
line once. This means that the amplitude can be written 
as 



W 



n A ^(^nj 



(65) 



where W is a set of c-webs that partition the 7i-particle 
state, and w is an individual c-web with n w particles. 
N\y is a combinatoric factor from exponentiation of the 
matrix elements and A w is the matrix element for the 
c-web w. 

This factorization is useful because collinear singulari- 
ties between final-state partons only exist when they are 
in the same c-web. Therefore if two partons in different 
c-webs are clustered into the same jet then the rate is sup- 
pressed by the area of the jet in rapidity and azimuthal 
angle, which is 0(R 2 ). At leading power, this implies 
that each jet's constituents are in the same c-web, and 
for a given set of c-webs W the measurement function 
can be factorized over this set: 



(66) 



Note that M. nm is the same measurement function as M. n 
but over <&„ iu instead of $„. Using the exponentiation of 
matrix elements, 



5(p^ ut ) = exp 



c-webs w 



where 



S w = J d$„„ A{$ n JM nw ($„„, pT 



(67) 



(68) 



As argued in Appendix |Aj RG invariance implies that 
the clustering logarithms can enter into at most the non- 
cusp anomalous dimension for the beam and soft func- 
tions at each order. For the px jet veto, the clustering 
logarithms will be a part of the v non-cusp anomalous 
dimension. 

Just as they did at 0(a 2 ), collinear divergences be- 
tween particles in the c-web will produce logarithms of R 
as a finite artifact of the divergence. Since there are at 
most n w — 1 collinear divergences in a c-web of n w par- 
ticles, this implies there are at most n w — 1 logarithms 
of R. Therefore, the general form of the clustering loga- 
rithms in the soft function S w for the c-web is, in the T 
veto case, 



C^ ra) (lni?) In 



A^(T cut ) = 
and in the pt veto case, 

AS w (pr) = (-Y W C&»>(lniZ) In 



/' 



(69) 



(70) 
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FIG. 3: Numerical size of various 0{a 2 s ) terms for a cut on prj (left panel) and on Tj (right panel) for vnu = 125 GeV at a 
8 TeV LHC. The terms are scaled by the leading order cross section, olo- We show two values of the jet radius, R = 0.4 (solid 
curves) and 7? = 1.0 (dashed curves) for the &t class of jet algorithms. Shown are the 0(a 2 ) contribution to the total inclusive 
cross section with no veto, a (2) (oo), the full 0(a 2 ) contribution with a jet veto, cr (2) (fc cut ) with fc cut = p™* or T cut , the 0{a 2 ) 
clustering logarithm corrections, Act' 2 ' (In R) , and the soft-collinear mixing contribution, Sa^(R 2 ). At small R, the clustering 
logarithms are sizable while the soft-collinear mixing terms are small. At large R the situation is reversed. 



Both C^p^ and Cp™ w ^ contain at most n w — 1 powers 
of lni?. These give rise to non-cusp anomalous dimen- 
sion contributions. RG evolution of the soft function and 
beam function clustering logarithms from the c-web w 
gives rise to an all-orders contribution to a of the form 



C / a™ 1 " ) (ms,Mj) = exp 
U^ w \vs,vj) = exp 





dfi 


( a s (fi) 


'Us 


A* 




j-Vj 


Av i 




1 Us 


v \ 





C^\hxR) 



R) . 
(71) 



bmce contains a In"™ 1 R term, when R ~ A 

this series is NLL. These terms are not directly obtained 
by resummation of lower order terms, underscoring the 
breakdown of resummation unless the coefficients of the 
clustering logarithms can all be calculated at once. 



V. CONCLUSIONS 

Jet vetoes are an important part of many experimental 
analyses. To obtain precise theoretical predictions, the 
resummation of large logarithms of the jet-veto scale is 
important, but they must account for effects of the jet 
algorithm. The algorithm generates two effects that can 
inhibit resummation, namely soft-collinear mixing and 
clustering logarithms. From a formal perspective, there 
are two scaling limits to consider, R ~ A and R ^S> A. 
Soft-collinear mixing arises at large values of the jet ra- 
dius R 3> A from correlations between soft and collinear 
radiation into the same jet. This mixing affects the re- 
summation of logarithms of the veto scale at NNLL and 
beyond and must be accounted for in order to prove an 



all-orders factorization theorem. When R ~ A, soft- 
collinear mixing effects are power suppressed, but clus- 
tering logarithms are important, which are a remnant 
of collinear singularities between particles in the jet that 
cancel in the total cross section. In the small R limit, new 
clustering logarithms at NLL arise at each order and can- 
not currently be resummed. Therefore, there is a tradeoff 
between these effects depending on the value of R cho- 
sen, and it is important to understand the relative sizes 
of these effects and incorporate them into estimates of 
theory uncertainties. 

In Sec. |TTJ we contrasted the factorization picture for 
inclusive observables Tb and Et with the exclusive jet- 
based observables 7} and pxj ■ The factorization and re- 
summation for the inclusive observables are well under- 
stood. In this paper we have explored the effects of soft- 
collinear mixing and clustering on the factorization and 
resummation for the exclusive observables. In Sees. Mil 
and |IV| we calculated the soft-collinear mixing and clus- 
tering logarithms at 0(a 2 ). 

The resummation for a pxj veto was recently consid- 
ered in Refs. []3J[T5]. In Ref. [TJ], the resummation was 
performed to NLL, but the clustering logarithms were 
only included at fixed 0(a 2 ). Our calculation of the 
leading clustering logarithms agree with the analytic re- 
sults given in Ref. [H] . Reference [TS] works in the limit 
A -C R ~ 1- They argue that in this limit the factor- 
ization formula in Eq. (32 1 holds up to 0(A) power cor- 
rections, and based on this extend the resummation to 
NNLL. While clustering logarithms do not pose a prob- 
lem for R ~ 1, we disagree that the factorization formula 
in Eq. p2| holds in this limit. Since soft-collinear mix- 
ing contributions are in fact not power suppressed for 
R ~ 1, the factorization formula does not reproduce the 
full NNLL structure of the cross section at leading power. 
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Therefore, the impact of these terms on the perturbative 
uncertainties found in Ref. |15j should be examined. 

We can investigate the numerical size of the jet- 
algorithm effects using Higgs production through gluon 
fusion as an example. The cross section for a veto k cut 
on Tj or pxj can be expanded in fixed order as 

a{k mt ) = <j lo + a« (fc cut ) + ct< 2 ) (fc cut ) + 0(a 3 s ) , (72) 

where o-M is C(a") relative to olo- The 0(a s ) terms, 
cr*- 1 -*, do not depend on the jet algorithm, and hence are 
the same for Tj,prj and Tb, Et, respectively. The 0(a 2 ) 
terms, d*- 2 ', include the effect of clustering logarithms, 
Aer( 2 )(lni?), and soft-collinear mixing, 5<t( 2 \R 2 ), 

a (2) (fc cut ) D Aa^QnR) , 5a^{R 2 ) . (73) 

In Fig. [3j we plot these terms in units of ctlo as a 
function of p™' (left panel) and 7~ cut (right panel), for 
R = 0.4 (solid curves) and R — 1.0 (dashed curves). 
For comparison, we also show the full 0(a 2 ) correction 
to the vetoed cross section, a^(k cut ), as well as the 
0(a 2 s ) correction to the total cross section without a veto, 

a (2)( fc cut = g7J|52]. 

The tradeoff between soft-collinear mixing and clus- 
tering logarithms is clear in Fig. [3] For R = 0.4, the 
clustering logarithms dominate over the mixing terms. 
For R = 1.0, the mixing terms dominate over the clus- 
tering logarithms. In each case, the size of the numeri- 
cally more important term is an appreciable fraction of 
the 0(a 2 ) correction. The difference between a^(k cut ) 
and a( 2 '(oo) primarily comes from the logarithmic terms, 
while the absolute size of <jW{k cut ) relative to <j( 2 )(oo) 
is indicative of the substantial cancellation between the 
jet-veto logarithms and the large NNLO if-factor for 
the total cross section 53J. Given that depending on 
R the clustering or soft-collinear mixing terms are a siz- 
able fraction of the logarithms ln(mir//c cut ) that one is 
trying to resum, it is crucial to understand their size at 
higher orders and correctly fold them into resummed un- 
certainty estimates. This has not been done in previous 
studies of the Higgs cross section where resummation for 
the pt veto is performed. 

The standard jet radii in Higgs analyses are R — 0.4 
for ATLAS and R = 0.5 for CMS. As seen in Fig. |3j 
with these values, at 0(a 2 ) the clustering logarithms 
are numerically important and the mixing terms can be 
regarded as power corrections. This suggests that the 
phenomenologically relevant limit is R ~ A, where soft- 
collinear factorization can be applied but clustering loga- 
rithms should be regarded as NLL. To proceed one should 
study the impact of the clustering logarithms on the re- 
summed perturbative series in order to properly take 
them into account in a resummed uncertainty estimate. 
For this purpose, a calculation of the O(ctf) clustering 
logarithms would provide very useful information. 

Note added in proofs: After this work was completed, 
Ref. |54j appeared, where a NNLL resummation formula 
for the cross section with a px jet veto is presented for 



R ~ 1 that accounts for the contributions due to soft- 
collinear mixing. This formula is equivalent to modifying 
Eq. ( 32 ) to include these terms by hand, and is consistent 



with our analysis that Eq. ( 32 ) by itself is insufficient 
to perform the resummation beyond NLL when R ~ 1. 
(Ref. [54] does not discuss the case R ~ A.) The structure 
of the 0(a 2 ) soft-collinear mixing terms suggests that for 
the pt case one may be able to absorb them into the beam 



and soft functions in Eq. ( 32 ) by performing additional 
collinear zero-bin subtractions in the soft sector. It is an 
open question and beyond the scope of this work if such 
an approach can be extended to all orders in a s such that 
Eq. ( 32 ) becomes a valid all-order factorization theorem. 
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Appendix A: Renormalization Group Constraints 

The anomalous dimensions of the beam and soft func- 
tions for the four observables that we study are con- 
strained by RG invariance. We give a brief summary 
of the properties of the anomalous dimensions and the 
general form for each observable. For the exclusive ob- 
servables, we work in terms of the veto variable 7~ cut or 
p™*, and these forms are useful when discussing cluster- 



ing logarithms in Sec. IV 



The hard function is universal for all the observables 
used, and the hard anomalous dimension is 



1h ( m H , M) = 2rs usp [a s (fj,)] In 2 



H 



2 7lrK(M)] ■ (Al) 



Consistency of the factorization theorem requires that 
the hard, beam, and soft anomalous dimensions cancel. 
For the four observables, this consistency is slightly dif- 
ferent. 

For beam thrust, the consistency relation is 

jff(m H ,^)S(TB) + 2m H ~f 9 B (m H TB 1 Li) + 7s'(7B,/ i ) = 0. 

(A2) 

Analogously, for the Tj < 7~ cut veto on jets, the consis- 
tency relation is 

7 f(m ff , M ) + 2 7 ^(m ff r cut ,M)+7f (T cut , M ) =0. (A3) 
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Both of these factorization theorems are described by 
SCETi. The scale dependence in the beam and soft func- 
tions is fixed by power counting; the beam function scale 
is fx 2 B = mjiT and the soft function scale is jxs = T . 
where T = Tb or T cut . This constrains the beam and 
soft anomalous dimensions, so that the coefficients of the 
cusp anomalous dimension are fixed relative to that of 
the hard function. For beam thrust, these anomalous 
dimensions are 

7*(t,M) = -2r» usp [a s (M)] 4a>(4) +7b[«.(m)] <*(*), 

7 f (k, fx) - 4r? usp [a s ( M )] U (£) + 7 f MA*)] *(*) . 

(A4) 

where £o(#) = + denotes the usual plus distri- 

bution. For a 7j veto on jets, the anomalous dimensions 
are 

7^(^r cut , M ) = -2r? usp [a s ( M )] In -JLL- 
+ 7sjctK(M)] . 

'T'CUt 

7 f(r cut ,M)-4r? usp [a s ( M )]ln^- 

+ 7l1etK(i»)]. ' (A5) 

The non-cusp anomalous dimensions for these two ob- 
servables agree at 0(a s ), but due to the different struc- 
ture of their factorization theorems (beam thrust con- 
tains a convolution while the Tj veto does not), the re- 
summed distributions only coincide at leading logarith- 
mic order. 

For the pr-based observables, which are described by 
SCETn, rapidity divergences exist and must be sepa- 
rately regulated. The rapidity renormalization group for- 
malizes the regularization of rapidity divergences by in- 
troducing separate anomalous dimensions for the beam 
(or jet) and soft functions [301 EH- This rapidity renor- 
malization group works much like the traditional RGE, 
with a scale v that is the analog to the usual /i scale. 

Power counting and renormalization group invariance 
provide strong constraints on the anomalous dimensions 
for py-based observables. Rapidity divergences are reg- 
ulated by factors inserted into soft and collinear Wilson 
lines; in each sector the regulator and its scaling are 

beam: z/"|n • P g \^ ~ i/'(A°)~'' , 
soft: ^\P 3g r' ^^(\y\ (A6) 

where P M is the momentum operator, and the group mo- 
mentum (P g ) of connected webs of gluons is regulated. 
In each beam function the large component of momen- 
tum is regulated. This power counting implies that in 
the rapidity RG the beam function scale is Vb = tuh 
and the soft function scale z/g — Pt, where pt = Et or 
p™'. In addition, in the standard RG in /x, the beam and 
soft functions live at the low scale, lib — Ms = Pt- 



For thepT-based observables, there are two consistency 
relations: one in /z-space and one in v-sp&ce. For the 
inclusive Et observable, the consistency relations are 

= 7» fl (m H , (x)5(E T ) + 2^ 9 (m H ,E T , li, v) 

+ 1 % 9B {E T , l x,v), 
= 2Y B 9 (E T , li, v) + i"s"{E T , li, v ) . (A7) 

And for a pTj veto on jets, the consistency relations are 

- -y 9 H 9 (m H , M ) + 2^ 9 {m H ,p c T u \ fx, v) 

= 2Yb 9 (Pt\ v) + i7 9 (Pt\^ v) ■ (A8) 

The general forms of the beam and soft anomalous di- 
mensions can be constrained through these relations as 
well as the exactness condition, 

^M = ^M, (A9) 

for F — B, S. We ignore running coupling effects in this 
exactness relation that will contribute to the v anomalous 
dimensions starting at 0(a 2 s ). For Et, the fx anomalous 
dimensions have the form 

j B ' g (m H ,E T ,fx,i>) 

= Ur^a^fx)} In — + 7 £> S (A*)] Wt) 

+ 7 > s ( M )]-£ (— ), 

// \ fX / 

= {^[^WJln^ + ^ 9 [a s {fx)]^5{E T ) 
-2^[a s (fi)]-cJ E ^ ), (A10) 

fX V LI J 

and the v anomalous dimensions have the form 

1b *{E t ,li,v) = -2T 9 usp [a s (fi)} ^o(^) 
+ -f*[a,(j i )]8(E T ), 

1 " s < a (E T ,(x,v)=4rZ uap [a s (n)] ^o(^) 

-Zf*[a a (n)]5(Er). (All) 

The cusp anomalous dimension dependence in the hard 
function fixes every part of the beam and soft anoma- 
lous dimensions except for the non-cusp terms and a cusp 
term 7 M [a s (fx)] in the beam and soft fx anomalous dimen- 
sions that cancels between them. Because the beam and 
soft functions have the same fx scale, as long as they are 
RG evolved in fx in the same way this 7 ^ term will can- 
cel. At fixed-order, this cancellation is guaranteed and 
has no effect on the fixed-order logarithms. It is ab- 
sent in the anomalous dimensions at 0(a s ), and may be 



15 



absent at all orders. Additionally, there is a constraint 
on the non-cusp anomalous dimensions from Eq. ( A7 1 , 



2 7fr[a s ] +2lB 9 l a s] + 7s "i a s} = 0. If the beam and soft 
functions are RG evolved in /i in the same way, then it is 
irrelevant how the \i non-cusp anomalous dimensions are 
divided between them, since only the sum matters. 

For the pxj veto the /j, anomalous dimensions have the 
form 



7£?etM^)]+7 AI MM)]m 



cut 
Px 



tS'W, v) = ^ cusp [a s ^)} In ^ + <y$f°[a a (n) 



2 7 > s (^)]ln 



cut 
Px 



(A12) 



and the v anomalous dimensions are 



cut 

Pt 



cut 

, M , „) = 4T| usp [a s (,i)] In ^ - Zy™[a.Qi)] . 

fx 

(A13) 

The same considerations as for the Ex observable apply 
here. All the cusp parts arc fixed by RG invariance with 
the hard function, and the only nontrivial unconstrained 
terms are the non-cusp v anomalous dimensions of the 
beam and soft functions. 

If clustering effects are associated with divergences in 
the beam and soft functions, then they can only impact 
the non-cusp beam and soft anomalous dimensions. This 
implies that clustering corrections have at most a single 
logarithm of the veto variable at each order. However, as 
we saw in Sec. |IV[ they can produce additional logarithms 
of the jet radius R, because these are not associated with 
any scale and are therefore not constrained by the RGE. 



Appendix B: Soft-Collinear Mixing at 0(a1) 



as A for them to be formally power suppressed. We also 
calculate the 0(a 2 s ) clustering contribution from inde- 
pendent emissions in the soft function, Aer^ cp ^ . 

The bare soft functions at 0(a s ) and 0(a 2 ) for inde- 
pendent emission in MS (without a measurement func- 
tion) are 



h{1) = as CA e^^ 2e r dfc+dfc - (fc + fc -)-i-. 

7T r(l - e) J V ; 



lr -b(i)ia 



(B2) 



where the label b stands for bare. We have implemented 
the on-shell conditions and left only the +, — light-cone 
components unintegrated. 

The gluon beam function B g can be written as the con- 
volution between a perturbative function I(t,z,(x) and 
the parton distribution function [5| 128] . 

B g (t,x,n)= I — , (B3) 

3={g,q,q} 

where t = rriHk + is the spacelike virtuality of the gluon 
entering in the hard interaction. The tree level, B g °\ 
bare 0(a s ) naive, Bg ^ , and zero-bin, B^ 5 ' , beam func- 
tions (without a measurement function) are [5] 

1 



B (0) 



,„ S(l,--m\ z) 



B 



b(i) _ 



OsCa 

2tt 



9(z)P gg (z) 



1 e 7BC // 



m H r(l - e) J 



dq + dq 



B 



b(i) _ (x s Ca 1 



9,0 



6(1 -z) 



7r rriH r(i — e ) 

/>00 

/ dq+dq- (q+r^iq-)- 1 -* . 
Jo 



ob(l) = Sb(l) 

g 9 



B 



b(i) 
9.0 ' 



(B4) 



In this appendix we present the 0(a 2 s ) calculations of 
the soft-collinear mixing contribution to the cross section, 



denoted as So^Sq^ 



- Sa^^P\ As explained in 



Sec. Ill this contribution arises from a correlation in the 



measurement between single independent emissions from 
the soft and collinear sectors, and formally breaks factor- 
ization if it is not power suppressed. These mixing terms 
have the same form as the independent emission cluster- 
ing terms at 0(a 2 ) in the soft and collinear sectors, and 
the total 0(a 2 s ) clustering correction to the independent 
emission contribution to the cross section is given by [see 



Eq. (47)] 



A(T indcp (2) = A(7 mdcp (2) + ^macp w + ^ W _ ^ 

The crucial result is that the soft-collinear mixing terms 
are NNLL terms that scale as R 2 , requiring R to scale 



indep (2) 



Jot (2) 



where 



Pgg(z) 



1 - z 



1 - z 



+ z{l-z) 



(B5) 



We have dropped the quark contributions to the beam 
function as they will not produce logarithms in the mix- 
ing terms. In B g ^ , there is a singularity as z — > 1 in the 
unregularized splitting function P gg that is regulated by 
q~~ . When considering the full measurement (both the 
factorized and soft-collinear mixing terms), there is an 
IR divergence in I that is canceled by the parton dis- 
tribution functions; however, this cancellation occurs in 
the factorized measurement function and is not present 
in the mixing term. The soft-collinear mixing correc- 
tion has a single divergence and B g ^' is proportional to 
(5(1 — z), meaning there is no nontrivial convolution with 
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the parton distribution functions. Therefore we will just 
compute the coefficient of the leading order cross section, 
ctlo: which is given by 

o-LO = cr H]$(m,H,li) J dYf g (x a , (i) f g (x b , n) , (B6) 

with x a ,b = (m H /E cm )e ±Y . 



1. Clustering Effects for Independent Emissions in 
the Soft Function 



We perform the soft function calculations first, as the 
techniques and results will be used later in the soft- 
collinear mixing terms. The independent emission terms 
are of the form 



0(Ti + T 2 < T cut ) - 0(71 < T cut )9(T 2 < T cut ) 



Only the integral over y t produces a divergence, and the 
integrals are straightforward to evaluate. The divergent 
term is 



a indep (2) rq-cut \ _ „ 1 ( V / H \ 4e 7T 2 2 

(B12) 

For PTj, we need to regulate the rapidity divergences, for 
which we use the analytic regulator jUj . In this case, the 
regulator factor is 



u 2a {k+k+)~ a = v 2a { VTlVT2 y a e^ 



(B13) 



The amplitude and measurement function are indepen- 
dent of yt- Hence, the soft function is proportional to 



SS V ' Jo * / dy t e 4ay * = dy t e 4a yt + 

J — oo J — oo J X 



The integral over A0 does not need to be regulated since 
there is no collinear divergence between the final-state 
particles. The measurement function for a Tj veto is 



4a 



dy t e 4 ^ 
0. (B14) 



AMi et(2) (T cut ) = 26»(Ai? 12 < R)\e(k+ + k+ <T 



9{k+ < T cut )9(k+ < T c 



(B8) 



We have multiplied by 2 to account for the case where 
the two gluons are both in the other hemisphere (where 
T is equal to the minus component of momenta), and 
dropped the (D(R A ) region of phase space where the two 
gluons are in different hemispheres but still cluster. For 
a PTj veto, the measurement function is 

AM^^ip^) = 29(AR 12 < R)[0(pri +PT2 < P c T ut ) 
- 9( PT1 < P c T ut )9( PT2 < p c T ut )] . (B9) 

In this case we drop the 0(R 4 ) correction from using 
the vector sum over transverse momenta in the combined 
constraint. 

To evaluate the matrix element, we use the variables 



We made the reason why this integral is zero explicit by 
breaking the range of integration into two parts: (—00, x) 
and {x, 00). In the lower range a > regulates the inte- 
gral and in the upper range a < regulates the integral, 
and the two terms cancel. (This is precisely equivalent to 
the case of a scaleless integral in pure dimensional regu- 
larization with a here playing the role of e.) Thus, 



Act 



indep(2)^cut>, _ 
SS 



0. 



(B15) 



2. Soft-Collinear Mixing Terms 



The mixing term between soft and n a -collinear emis- 
sions at 0(a 2 ) is given by 



Jo n 



jet (2) 
as 1 



(B16) 





Ay = 




k x k 2 


yt 




fcj k 2 




k 1 k 2 


k 1 k 2 


-y-cut 


: Ti = 


k+ 




r 2 


- k+ 

— ^2 > 




Pt ■ 


PTl = 




~k~ 

"a j 


Pt 2 




k 2 . (BIO) 



For the Ti veto, the soft-function contribution is 



Aa mdcp(2) (rcut) 

vlo /« s C A \ 2 (e 7E ^ 2 ) 2e 



e(AR 12 < R) f°d7i dT 2 (Tr 75)- 
Jo n Jo 



(Bll) 
dAy 



and similarly for (SCb). We have included the integral 
over the relative azimuthal angle A(f> between the soft and 
collinear particles. Since we are concerned with the lead- 
ing divergences, we do not need to regulate <j>. The 0{a 2 ) 
measurement corrections for 7} and pTj, SAi^ ct< ^ 2 \T cut ) 
and SM^^d^), are 

5M' sct (2) (T cut ) = 9(AR SC < R) \9{T C + T< T cut ) 

-9(T c <T cut )9(T s <T cut )], 
5M ict (2) (Pr ut ) = 0(AR sc < R) [9(p Tc +PTs < Pt") 

-9(p Tc <p c T ut )9( PTs <p T " t ) . 

(B17) 
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For the Tj veto, the two soft-collinear mixing terms, 
(SC a ) and (SCb), are equal. They are made up of the 
naive and zero-bin beam function contributions, so that 



(SC) = (SC) - (5C (0 »). 



(B18) 



The zero-bin of the beam function in each collinear sector 
is equal to the soft function contribution, meaning 



(5C (0) ) = 2(55). 
Changing variables from k~ to Ay, 



(B19) 



q + z 



m H (l - z) 



the total correction from the naive beam function term 
in the mixing is 



(5C*)(T cut ) 

'a s C A \2(e 7 V) 2e 



7T J r(i-e) 2 /; " 
rdA</> 



z \ 2e 1 



1-W m^ +2e 



/ dAye 2 ^ / -_-0(Ai?<i?) 

J-oo JO ^* 



x ; dq + dk + (q + k^ 
'o 

•»+ ^ T"cut 



-l-2e 



(B21) 



g+ < r cut ) - %+ < r cut )6»(fc+ < r cut ) 



The remaining integrals are finite, so we can set e = 
in the integrand after pulling out an overall scale depen- 
dence of (T cut ) -2e . The integrals over Ay and A<fi give 



pdA„ r dA ^^<i 2) =ii? 2 

J-oo JO T 2 

and the integrals over q + and fc + give 
' [ 6{k+ + q+ < r cut ) 

+ < .put) < r cut^ 



o 



Thus, 

(5<5)(T cut ) 



it / m# VI — z 



(B22) 



(B23) 
"6 ' 

2 c 



V mfl-7" cut / 6 



2 « 2 



(B24) 



Expanding z as a distribution about 2 = 1 and keeping 
only the divergent term, 

(j(z)P S3 (z)( r ^-) 2£ = -U(\ -z)+ O(e ) , (B25) 

the soft-collinear mixing term is 



(B20) (5C)(T cut ) = 







( , 2 \ 






{m H T cut J 



2 ^2 



G 



-i? 2 



(B26) 

Subtracting the zero-bin contribution, which is twice 
Eq. (B12|, the correction to the cross section becomes 



(-4 



\m H T cut J \(T cut ) 2 
\ 7T J 3 T cut 



2 c 



(B27) 



For the pxj veto, the soft independent emission, and 
therefore the zero-bin contribution, are zero. Hence, the 
soft-collinear mixing terms are just given by the naive 
contribution. To evaluate the matrix element, we change 
variables from q + , k + to PtciPts, where 



PTc z 



H -—7^-,, k+= PTs e-y^^±. (B28) 
m H (l ~ z) p Tc z 



We also use the same change of variables from fc~ to Ay 
as in the Tj veto. There are again rapidity divergences 
not regulated by e for which we use the analytic regulator. 
The regulator factor for (SC a ) is 



v 2a {q + k+)- a = v la { V TcVTsY a 



PTcZ 



m H {l - z) 



(B29) 

The integrals over Ay and A<fi can be performed as be- 
fore, and we find 



C a\ 2 (e~< E u 2 ) 2e 1 - /l 



(1 - e) 2 m H 
dp Tc dpTs Pt, 



2a. 



l-2e-3ct -l-2e-a 
PTs 



(um H ) 2a R 2 



e( P Tc + PTs < Pt") - V(PTc < Pt") 0(pt s < Pt") 



(B30) 



We can set e = everywhere, as there is only a rapidity divergence as z — > 1. Performing the final integrals, we 
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find 



vrriH 
(p° ut ) 2 



2q ^2 



-1 - z\ 2 " 

z 



-R z 



(B31) 



Expanding in z to extract the divergence, we get 



i? 2 



(pip*) 2 J 6 

(B32) 

This gives the following correction to the cross section 



6a 



J ot ( 2 ) ( m CUt 

SC a 



i cutN 1 (a s C A \ 2 ( vm H \ 2a 7r 2 



(i 



i? 2 . 
(B33) 

The regulator for the mixing between soft and n&- 
collincar sectors gives a difference scale dependence. In 
this case the label component of collinear momentum is 
regulated, which effectively amounts to regulating the 
minus momentum component in the above calculation. 
This changes the scale dependence to (v/mH) 2a with an 
overall minus sign: 



■(v/m.ff) 



2<< 



(B34) 



Thus the entire soft-collinear mixing contribution to the 
cross section, 8cr 



jet 

sc a 



So £* is 



)2ot / \ 2a. 

- i - i 
\m H ) 

fa s C A \ 2 2n 2 2 m H 

-ctlo -=--« In-^ 

V 7T / 3 p™ 



(B35) 



i^T measurements. These measurements are IR safe, and 
when a set of particles becomes collinear their contri- 
bution to the primary measurement is the same as the 
contribution to the observable that is vetoed (the px or 
T of the jet). Here, we only keep the divergent term with 
the appropriate p™* and 7~ cu * scale dependence, as this is 
the piece connected to the soft and beam functions. The 
finite terms do not take part in the resummation and can 
be captured as usual by matching the resummed result 
to the full fixed-order result at NNLO. 

The 0(a 2 ) soft measurement function for the cluster- 
ing correction relative to the inclusive measurements is 
[see Eq. poll] 



AM{ cti2) {k cut ) 



9(AR 12 < R)0(k iet < k cut ) 

6(AR 12 > R)9(h < k cut )9(k 2 < fc cut ) 
8(ki +k 2 < fc cut ) , (CI) 



where k = T or px- In this measurement function, fcj e t 
is the observable for the clustered pair. For k = T, 
7jot = T\ + T 2 except when the jet spans the bound- 
ary at the Higgs rapidity. For k — px, the scalar sum of 
transverse momenta is not the same as the magnitude of 
the vector sum. However, when R ~ A, these differences 
become power-suppressed for both the T and px observ- 
ables. Since this is the limit we are working in, and we 
are neglecting power-suppressed terms in this limit, the 
measurement functions for the clustering effect simplify 
to 

AX jct(2) (fc cut ) = 6(AR 12 > R)\e(ki< k cut )6(k 2 < k cut ) 
- 9(h + k 2 < fc cut )l . (C2) 



Appendix C: Clustering Logarithms in the Soft 
Function at Ofa,) 

Clustering logarithms first occur at 0(a 2 ), and they 
are easiest to calculate in the soft function. RG invari- 
ance can be used to extract the beam function contribu- 
tions, as explained in Sec. |IV| 

The definition of what is a clustering effect is subtle 
because one must define what the effect is relative to. 
One ostensibly natural option is to define it relative to 
the cross section if no clustering takes place. However, as 



explained in Sec. IV this measurement is infrared unsafe 
due to collinear singularities between partons. A more 
sensible "primary" measurement are the inclusive Tb or 



The phase space constraint on Ai? 12 suggests a conve- 
nient set of coordinates for the calculation. In terms of 
the rapidity yi, azimuthal angle (f>i, and observable ki 
(where k = T or px), the coordinates we use are 



vt = 2(2/1 +V2) 
&y = yi - V2 , 
h — h + k 2 , 



1 



+ 61), 



A(f) = 4>i - <t>2 ) 

h 



(C3) 



The full 0(a 2 ) soft function matrix elements can be 
found in Ref. |55j . Using k = px, in terms of these coor- 
dinates the non-Abelian soft matrix elements arc 
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1 



1 



1 



A p^ t z 2 (l — z) 2 cosh Ay — cos A0 z 2 + (1 — z) 2 + 2z(l — z) cosh Ay 
x j(z 2 + (1 - z) 2 ) cos A<p - z(l - z)(l - cos Ac/) cosh Ay) 



z 2 (l - z) 2 sinh 2 Ay 



A f = Ag 4 C A T R n f 



(cosh Ay - cosA<^)(z 2 + (1 - z) 2 + 2z(l - z) cosh Ay) J ' 
1 1 1 



;4 t z 2 (l-z) 2 (cosh Ay- cos A0) (z 2 + (1 - z) 2 + 2z(l - z) cosh Ay) 



x <^ z(l - z) - 



2z 2 (l - z) 2 sinh 2 Ay 



(cosh Ay - cos A</>)(z 2 + (1 - z) 2 + 2z(l - z) cosh Ay) 
In the small AR limit, the matrix elements simplify to 

1 2 



(C4) 



A% = 4g*C A 



Af = 4g A C A T R n f 



p% t z 2 (\-z) 2 AR 2 
1 



2 ,i 2z 2 (l-z) 2 Ay 2 

- + ^-^ + AR 2 

4z 2 (l-z) 2 Ay 2 



z(l-z)- 



AR 2 



p% t z 2 {\- z) 2 AR 2 

We start with the calculation for pTj ■ The measurement function can be written as 

1 



AMf'ip^) =0(AR> R)0 



p C T Ut < Pit < FT' 



max(z, 1 — z) 



(C5) 



(C6) 



The matrix element and measurement function are independent of y t , meaning the integral over y t is unregulated. 
To regulate this rapidity divergence we use the rapidity regulator, which regulates the z-component of the group 
momentum for the c-web through the factor 



Prt 



2zsinh(y 4 + ^ Ay) +2(1 - z)sinh(y t - *Ay) 



Integrating over y t then gives a single l/r] divergence and this is the only divergent part of the calculation, 
i/^I' J dy t |2zsinh(y t + ^Ay) + 2(1 - z) sinh(y t - ^Ay) |"" = Vf>"? + 0{if ) . 



(C7) 



(C8) 



The matrix element scales simply with prt since it is the only dimensionful variable, and so prt can be easily integrated 
against the measurement function including the regulator factor. This integral is 



/>oo 

/ dpTtp"^'* 1 6 
Jo 



cut ^ , »,cut 

Pt < Prt < Pt 



-(PtT 4 "" 7 ' ln[max(z, 1 - z)] . 



max(z, 1 — z) 

Carrying through the integrals over the on-shell conditions, the soft function correction for the pt veto is 

n ,i ,oo r n dA(j) ln [ max ( /2) ! _ z y 



(C9) 



z(l-z) 



- 6(AR > R) p\ t z 2 {\ - z) 2 A(z, Ay, A<j>) 



The remaining integrals are finite and have the form 

a\nR + b + 0(R) . 



(CIO) 



(Cll) 



We determine a analytically and extract b numerically. To determine the coefficient of In R, we rewrite the matrix 
element as the difference 



A = (A - A R ) + A 1 



(C12) 
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where A R is the matrix element expanded in the small R limit, given in Eq. (C5). The difference A — A R is finite as 
R 0, meaning that up to 0(R 2 ) corrections, 



JdAydAct>0(AR> R) (A-A R ) = J dAydA<j)(A-A R ) + 0{R 2 ) 



(C13) 



Integrating the matrix element A R is simple, using the relations 

r7! dAcf, 1 



dAy 

O 

dAy 



7T AR 2 
" dAcj) 2Ay 2 



6(AR >R) = -lni? + ln27r + C(i?) . 

f 



, „. 6(AR > R) = -l nj R + ln27r+ - +0(R) 
7T Art 2 



(C14) 



Carrying through the integrals, we obtain 



dz 



dAy 



dAtfi ln[max(z, 1 — z)] 



z(l-z) 



V 4 Tt z 2 {\- z) 2 A\ 6{AR> R) 

4 ff 4 Cl{-^(131 - 12tt 2 - 1321n2)lni?+ — [-13 + 121n2 + 21n(27r)(131 - 12tt 2 - 132 In 2)]} 
dA(/> ln[max(z, 1 — z)\ 



dAy 



z(l-z) 



V% t z 2 (\~ z) 2 A R 9(AR>R) 



4:g 4 C A T R n f ( — (23 - 24 In 2) In R + — [13 - 12 In 2 - 2 1h(2tt) (23 - 24 In 2)1 } 
I- 18 36 J 



18 v ' 36 

Performing the full integrals numerically, we get 



(C15) 



dz 



dAy 
dAy 



*dA(j) ln[max(z,l - z)] 



7T z(l — z) 

* dA<t> ln[max(z, 1 - z)] 



V % t z 2 {\ - z) 2 {A A - A%)\ = Ag A C 2 A (4.66) , 
p 4 Tt z 2 (l - z) 2 (A f - A?)] = Ag 4 C A T R n f (0.138) . 



7T z(l — Z) 

Thus, the final result for the leading divergent corrections due to clustering in the bare soft function is 



(C16) 



b (2) /nuts 



AS b(2) (p. 



1] V 7T 



( 7T ) V P ™* j { 



(131 - 12tt 2 - 132 In 2) InR - 1.12 



C A T R n f -^(23 - 241n2)lni? + 0.764 | 



(C17) 



Note that the constant terms, —1.12 and +0.764, depend on the choice of inclusive observable that the clustering 
effect is defined relative to. Hence, they are different from those found in Ref. [H], since here we use = |pri| + \Pti\ 
while Ref. [TJ] uses the vector sum |pVi + Vti \ ■ When using the latter in our calculation we reproduce the results in 
Ref. [H]. 

The steps to calculate the clustering logarithms for 7} proceed analogously. The calculation is slightly more tedious 
due to the fact that a particle's contribution to the observable changes depending on what hemisphere it is in. 
Therefore, there are two regions of rapidities j/i^ of the two particles to consider: 

i) 2/i j 2/2 > and 2/1,2/2 < 0: Both particles are in the same hemisphere. These configurations will contribute to 
the divergent terms. 

ii) yi > 0, j/2 < and y\ < 0, y% > 0: The particles are in opposite hemispheres. This region does not contribute to 
the divergent terms. The region of phase space where particles can cluster scales as 0(R 4 ) and can be neglected. 



The variables in Eq. (C3) are also useful for this calculation. Carrying through the calculation, we find the divergent 

){' 



terms contributing to the bare soft function 

A5 b(2) (T cut ) = l/a,\ 2 
4e V 7r / 



C A T R nj 



n~ cut 

- 1 

-~9 ( 
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'131 - 12tt 2 - 132 In 2) In R - 0.937 



23 - 24 In 2) In R + 0.747 



} 



(C18) 
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Note that the divergent mi? terms are the same as for prj, while the constant terms differ. 
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